FIG. 1 

Example Spectral Density Function 




frequency 



FIG. 3 

Example Amplitude Function Made from 10 Piecewise Continuous Functions 

(hollow dots are the knot points) 
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FIG. 4 

Example Amplitude Function Composed of Multiple Functions with Overlapping 

Domains 

(the top curve is the sum of the curves beneath it) 
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Figure 5 

Example Phase Angle with Doppler 




FIG. 6 

Example Phase Angle with Doppler when Emitted Frequency Selected to Produce 
Approximately Linear Phase Angle Response 



2.08P-K36 



frequency (Hz) 
7.dS^eM3B 2.084e-K)6 2.086.e40S 



2.08qB-K]6 2.09P-K36 2.092.e-K]6 



FIG. 7 

Example Phase Angle with Doppler when Emitted Frequency Selected to Produce 
Always Positive and Approximately Linear Phase Angle Response 

16-- 
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FIG. 8 

Example sin(Phase Angle) with Doppler when Emitted Frequency Selected to 
Produce Positive Phase Angle Response Close to Zero Phase Angle 
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FIG. 9 

Example cos(Phase Angle) with Doppler when Emitted Frequency Selected to 
Produce Positive Phase Angle Response Close to Zero Phase Angle 
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FIG. 10 

Example Phase Angle with Doppler using 2 MHz Emitted Signal 




FIG. 11 

Example cos(phase angle) with Doppler using 2 MHz Emitted Signal 
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FIG. 12 

Example sin(phase angle) with Doppler with 2 MHz Emitted Signal 




FIG. 13 

Six Parameter Distribution for Velocity 




velocity (m/sec) 
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FIG. 15 
Parameters for Example 



Number of function segments: N = 7 


M = 20 


Lowest frequency in basis: = 0.1996961309 


10^ Starting 




time: = 0.00006843506494 


Frequency increment: = 2207.908286 


Ending time: = 0.00007493194556 
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Figure 16 
Example ||A|| Matrix 
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FIG. 18 
Example Envelope Function 



Gaussian Pulse Envelope Function 
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FIG. 19 

Spectral Density Component Functions Used to Calculate Example f(t) Values 



Example Actual Spectral Density and Spectral Density Component Functions 
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Actual F(^ 
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FIG. 20 

Example Signal and Locations of Measurement Points 



g(t) and Points Where Measurements Taken 




FIG. 21 

Example Signal after Adjust for Gaussian Envelope and Locations of Measurement Points 
(The period of this signal is much longer than the short segment shown.) 
f(t) and Points Where Measurements Taken 




FIG. 22 

Example Calculated Results for G(f) and Actual 
G(f) 

(The calculated result overlays the actual G(f).) 

G(Q: Actual and Predicted Spectral Densities 




frequency 

Actual G(f) 

G(t): distribution calc'd usina l|PI| ■ 

o o o o o o Segment endpoints for precfictea ( 



FIG. 23 

Example Calculated Results for F(f) and Actual F(f) 
(The calculated result overlays the actual F(f).) 

F(f): Actual and Predicted Spectral Densities 




frequency 



F(f): distrib ution calc'd using ||P|| 

o o o o o o Segment endpoints for predicted 
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FIG. 26 

Three Frequency Function with Four Digitization Bins 



Actual and digitized results forf(t) = 2.5*sin(200*Pi*t)+10*sin(240*Pi*t)+4*sin(1050*Pi*t) 
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FIG. 27 

Three Frequency Function with Four Digitization Bins and Least Squared Errors Estimate 
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FIG. 28 



Three Frequency Function with Four Digitization Bins and 
Least Absolute Value Errors Estimate 



m 



- 


ej^ o 

r \ A 

0 p ooo o 


o 

A 

uooWo do oo 


oo 

P\o4 qOOOO oo 


• oldo5 


/vQ.bl I 0.015 


0.02^ \ 


0.p25 0.03 1 






time I y 




r ° 


t oo oVo<^o 


o P o/> 


So do b o^ ojo 








„ V 










J oo o* 




o 


0 o 



+ + * + f(t) -- before digitization 
o o o o o o f(t) -- after digitization 
rft) - before digitization 

Estimated f(t) using least absolute value errors after digitiz 



FIG. 29 

Three Frequency Function with Four Digitization Bins and 
Least Squared Errors and Least Absolute Value Errors Estimates 




f(t) » before digitization 
ffl] after digitization 
m - before digitization 

Estimated fO) using least squared errors after digitization 
Estimated f$[) using least absolute value erors alter digitizatioi 



FIG. 30 

Actual and Estimated G(f) Calculated without 
Accounting for Digitization 




Actual G(p 

G(f): distribution calc'd using ||P|| without grs 
Segnnent endpoints for predicted quadratic 



FIG. 31 

Actual and Estimated F(f) Calculated without Accounting 
for Digitization 

F(0: Actual and Predicted Spectral Densities 
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FIG. 37 

Double Sided Constraint Violation Value Function 
(Digitization bin width = 1 ) 
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FIG. 40 

Tapered Heaviside Gate Function Examples 

(k«1/2.p = 4| 




[k = 199/200. p = 4010] 
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FIG. 41 

Tapered Heaviside Gate Function for Extreme Value of G(f) is Less than 

Zero 

Function for sign of G(jl(fext). [k = 1/2. p = 8] 





FIG. 42 

Example G(f) Calculated Using NCVl Nonnegativity Constraint Value 

Function 

G(f): Actual and Predicted Spectral Densities 
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FIG. 44 Basic Calculation Sequence 
Page 1 of 3 



1. Specify solution parameters 


D. 


Select type of nonnegativity 


A. 


Specify measurements 




constraint (refer to FIG. 38 and 




a. Number of measurements 




FIG. 39) 




(M) 




a. Scaled by G/fj,ext)' 




b. Time interval between 




b. Scaled by -Gj(f;,ext) Select 




measurements 




peak count constiraint (refer 


B. 


Specify basis frequencies 




to FIG. 43) 




a. Lower and upper frequency 




a. Select whether will use 




limits (fL, fw) 




peak count constraints 




b. Number of frequency 




b. Select number of peaks for 




intervals (N) 




G(f), F(0 


C. 


Select Objective Function 


F. 


Specify properties of 




components (refer to FIGS. 33, 




digitization bin constraints 




34, 35) 




a. Sharpness of constraint 




a. Arc length measure for 




penalty (final value of 




G(f),F(f) 




parameter w) 




b. Bounded area measure for 








G(f), F(f) 




final calculations (final 




c. Quadratic curvature for 




value of parameter p) 




G(f),F(f) 


G. 


Specify rate of change of 








nonnegativity constraints 








(parameter k) 








proceed to precalculations 



FIG. 45 Basic Calculation Sequence 
Page 2 of 3 



2. Precalculations (before search for 


C. 


Calculate OF using starting |d| 


solution) 




and calculate starting OF scale 


A. 


Calculate matrix ||A|| (Refer to 




factor X (A=10^(order of 




FIG. 25) 




magnitude of OF minus 2) 


B. 


Calculate objective function 


D. 


Select starting value of 




(OF) expression and its 




parameter p, e.g., p=2. 




Hessian (this will be a matrix 


E. 


Calculate starting values for 




of numbers) 




constraint value functions: 


C. 


Calculate nonnegativity 




digitization bin constraint 




constraint value (NCV) 




violation value function 




function and its gradient and 




(CVV), nonnegativity 




Hessian (these will be a vector 




constraint value function 




and a matrix of expressions in 




(NCV), peak count constraint 




terms of |a| and |x|). 




value function (PCV) 


D. 


Calculate peak count constraint 


F. 


Calculate value for constrained 




value (PCV) function and its 




objective function (COF): 




gradient and Hessian (these 




COF= OF/X + CVV + NCV + 




will be a vector and a matrix of 




PCV 




expressions in terms of |a| and 








|x|). 




proceed to multidimensional 


3. Initialize 


Newton method search. 


A. 


Select starting |a| and |x| - e.g., 








constant vectors with small 








values 






B. 


Create starting |d| by stacking 








|a| over |x| 







FIG. 46 Basic Calculation Sequence 
Page 3 of 3 





4. Multidimensional Newton Method 


(4) IfCOFk+,<COFk then goto 




Search (designate starting value of 


step c, otherwise got to 


5 


variables with subscript k, e.g., 


step b.(2) and reduce |s|. 




starting |d| as |d,,|.) 


unless |s| has been reduced 




A. Using current value of p 


to a very small value, in 




a. (This sequence calculates 


which case goto to step d. 




new |d|) 


c. (Get here when have lower 


10 


Using |dj 


COF) If significant change in 




(1) Calculate gradient for 


|d,^+,| from |dk| then 




OF 


(1) set|d,Hd,,i| 




(2) Calculate gradient and 


(2) go back to step 4.A.a 




Hessian for CVV. 


otherwise go to step d. 


15 


(3) Calculate gradient and 


d. (Get here when stopped getting 




Hessian for NCV. 


changes in Id^+il) 




(4) Calculate gradient and 


If at final value of parameter p 




Hessian for PCV. 


then go to step B, otherwise 




(5) Sum gradients and sum 


(1) Increase p, e.g., 


20 


Hessians and use these 


multiply current value 




to calculate the change 


ofpby 3/2 




in|d|: |s|. 


(2) Calculate OF using \d^\ 




(6) Calculate |dk+,|=|dj + |s| 


and the new p 




(7) Calculate new OF, 


(3) Calculate new value 


25 


CVV, NCV, PCV and 


for scale factor X. 




. COF,,j. 


(4) Go back to step 4.A. 




b. (Backstep, if necessary) 


B. (Get here when have finished 




(1) IfCOFk^,<COFkthen 


calculations with final value of p) 




goto step c. 




30 


(2) Reduce |s| by order of 


** Search complete ** 




magnitude 






(3) Calculate |dk+i|=|dk| + |s| 





If 



1 



FIG. 47 

Actual and Calculated G(f) from Digitized Data 
M=50 

G(0: Actual and Predicted Spectral OensHlBS 
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FIG. 48 

Actual and Calculated F(f) from Digitized Data 
M=50 

F(0: Actual and Predicted Spectral Densities 
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FIG. 49 

Actual and Calculated Spectral Density Function from 
Digitized Data 
,1VI=50 
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FIG. 50 

Actual and Calculated G(f) from Digitized Data 
M=200 

G(f); Actual and Predicted Spectral Densities 
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FIG. 51 

Actual and Calculated F(f) from Digitized Data 
M=200 

F(f): Actual and Predicted Spectral Densities 
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FIG.52 

Actual and Calculated Spectral Density Function from 
Digitized Data 
M=200 

Total Spectral Density Function: Actual and Predicted Spectral Densities 
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FIG. 53 

Actual and Calculated Cumulative Spectral Density 
Function from Digitized Data 
M=200 

Tctil Anplhjda Cumulitiv« Ofstribulion Bistdon i Sptotnl Dtntity Function Using 7 Sagmwils irith 300 Musummnts. 




2e406 2.002e'K)6 . 2.004e-t€6 2.Q06e406 2.0069406 



frequency 



Estimated Cumulative Total Amplitude 
Actual Cumulative Total Amplitude 



